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Abstract 

The quantum diffusion of a particle in an initially localized state on a cyclic 
lattice with N sites is studied. Diffusion and reconstruction time are calcu- 
lated. Strong differences are found for even or odd number of sites and the 
limit N — > oo is studied. The predictions of the model could be tested with 
micro - and nanotechnology devices. 

I. INTRODUCTION 

The problem of a classical particle performing a random walk in various geometrical 
spaces 1 and the quantum random walk 2 have been thoroughly studied and compared. The 
classical and quantum cases have striking differences. One of these differences is that, 
whereas the classical spread increases with time as a/T, in the quantum case we have a 
stronger linear time dependence of the width of the distribution. In the quantum mechanical 
case, we can identify two different causes for the spreading of the probability distribution 
describing the position of a particle. There is a spreading of the distribution caused by 
the random walk itself, also present in the classical case, and superposed to it, there is 
the quantum mechanical spreading of the probability distribution due to the time evolution 
of a particle in a localized state. This second type of spreading is the main interest of 
this contribution. For this study, we will consider a quantum mechanical particle initially 
localized in one site of a one dimensional cyclic lattice with N points. In most treatments 
of quantum random walks in a lattice it is assumed that the number of sites, N, is large 
compared with the number of jumps of the time evolution and therefore the system does 
not notice whether the lattice is infinite or cyclic, that is, finite with periodic boundary 
conditions. In our analysis we will not assume that N is large and we will find some peculiar 
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and interesting features like, for instance, a very different behaviour for even or odd values 
of N. There are several motivations, besides the general academic interest, for allowing low 
values of N. For instance, cyclic lattices with a few sites have been built (quantum corrals) 
with nanofabrication techniques and in quantum computers we deal with systems with N = 2 
(qubits) or iV = 3 (qutrits). In all these cases we may be interested in the quantum diffusion 
time of an initially localized state. The quantum behaviour in the continuous case N — > oo 
is also interesting because it can be experimentally tested building small conducting rings 
with microfabrication techniques. 



II. DEFINITION OF THE MODEL 



In this work we will consider a particle moving in a one dimensional periodic lattice with 
N sites and lattice constant a represented in Figure 1. The quantum mechanical treatment 3 
of this system requires an N dimensional Hilbert space 7i. The lattice sites will be labelled 
by an index x running through the values 0, 1, • • • , N — 1. We will adopt a very useful 
notation for the principal N th root of the identity defined by 

uj = . (1) 

Integer powers of this quantity build a cyclic group with the important property 

1 =u Nn ,Vn = 0,±1,±2,-- - . (2) 

The position of the particle in the lattice can take any value (eigenvalue) a(x — j) where 
a has units of length, j = (N — l)/2 and the integer number x can take any value in the 
set {0, 1, • • • , N — 1}. The eigenvalues have been chosen in a way that position can take 
positive or negative values in the interval [—aj, aj] . Notice that j is integer for odd N and 
half-odd-integer if N is even. The state of the particle in each position is represented by a 
Hilbert space element tp x and the set {ip x } is a basis in 7i. In the spectral decomposition, 
we can write the position operator X as 

AT-l 

X = a(x - j)<p x (<p x , ■) , (3) 

x=0 

that clearly satisfies X(p x = a(x — j)tp x . Momentum is formalized in the Hilbert space by 
means of a basis {<f> p }, unbiased to the position basis, where p is an integer number that can 
take any value in the set {0, 1, • • • , N — 1}. The momentum operator is given in terms of its 
spectral decomposition as 
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N-1 

P= Y,9(p-j)M<i>pr) , (4) 

p=0 

where g is a constant with units of momentum. The eigenvalues of P have been defined 
in a way to allow for movement of the particle in both directions, anti-clockwise (positive 
eigenvalues) and clockwise (negative eigenvalues) along the circular lattice. Notice however 
that the state of zero momentum is only possible when N is odd. We will find in this 
work that there are several important differences in the system when N is even or odd. 
The position and momentum bases are related by a unitary transformation similar to the 
Discrete Fourier Transform 

1 N-1 



N p=0 
and 

-I N-1 



4>v = ~m E ub-M*-»+°l*-P) <p x , (6) 
ViV x=0 



{<p x , <f> p ) = -A^Jp-iKz-j) + «(z-P) > (7) 
v jV 



with 



where a is a parameter such that 

f for iV odd, 

a = 1 ( 8 ) 
I 2 for N even. 

The constants a and g are not independent but are related by 

agN = 2tt . (9) 

This condition follows 3 from the requirement that in the limit iV — > oo the commutation 
relation of position and momentum should be [X, P] — > % (we adopt units such that Ti — 1). 
Notice that for finite N, the commutation relation [X, P] — % is impossible. 

All these definitions are compatible with the physical requirement that momentum is the 
generator of translation and position generates increase in momentum. That is, 

e~ iaP ip x = u a ip [x+1] 

where the symbol [•] denotes modulo N, that is, [N] = 0. Notice that an iV-fold application 
of these translation operators is equal to the identity 1 if N is odd but is equal to — 1 if N 
is even. This is reminiscent of a 2n rotation of a spin 1/2 system. 
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In equations (5), (6) and (7) we could absorb the phases ax and ap in the bases {<p x } 
and {4> p } by an appropriate phase transformation (this is only relevant for even N because 
a^O). However this option would result in a complication of Eqs. (10) where the phase u a 
would not appear but a sign change would appear in the translation from site x = N — 1 to 
site x = and also from p = N — 1 to p = loosing thereby the homogeneity of the lattice 
because not all lattice sites would be equivalent. Later in this work it will be convenient to 
take this option. 



III. SPREADING OF A LOCALIZED STATE AND DIFFUSION TIME 

At any instant of time, the state of the particle will be determined by a Hilbert space 
element \&(t) that in the position representation is given by the coefficients c x (t) such that 

N-l 

*(*) = E c *(*K • (ii) 

A given state \l/(0) at an initial time t — will evolve according to the time evolution unitary 
operator given in terms of the hamiltonian H as 

U t = exp(-iHt) . (12) 

In this work we are interested in the time evolution of a state corresponding to a particle 
initially localized in a lattice site, say at x — 0, at rest, that is, with (P) = 0. Such a state 
is given by ^(0) = (po, that is, c x (0) = 5 Xj0 . Let us assume a free particle with hamiltonian 
H = P 2 j 1m. With this hamiltonian we can easily find that the state for any time will be 
given by 

i N-l 

C *( T ) = E J x{p - j+a) -^ 2T ) , (13) 
^ v p=0 

where we have introduced a dimensionless time parameter T = t/r with a time scale r 
defined by 

2ma . . 

r= — , (14) 

that, as we will see later, corresponds essentially to the diffusion time. We have chosen 
the free particle hamiltonian, however many of the following results do not depend on the 
specific form of this hamiltonian and are also valid for any hamiltonian invariant under the 
transformation P — > —P. During the time evolution of a particle, initially in the site x — 0, 
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the expectation value of the position and momentum will remain zero, (X) = (P) = but, 
due to the quantum spreading of the state, the probability distribution of the occupation 
of other lattice sites will grow. We will study some features of this quantum diffusion. The 
probability of occupation of the lattice site x at time T is given by 

\c x (T)\ 2 = \ YY, ' u (p ~ q)x - (p - q)(p+q - 2j)T . (15) 

One of the sums can be analytically performed after a change of the summation indices but 
it is not really convenient to do it. 

Due to the periodicity of the lattice, we expect that the amplitudes and probabilities 
of Eqs.(13) and (15) will be periodic in time. This is indeed the case but with different 
periodicity for N even or odd. That is, for the amplitude we have 

cAT) = l^ + N) Cor Note 
\c x (T + 4N) for N even 

and for the probability we get 

, cm ,2 = {MT + AOI 2 for TV odd 
" \ \c x (T + N/2)\ 2 for N even 

It is remarkable that the period of the amplitude is equal to the period of the probability 
for iV odd, but it is eight times longer if N is even. The periodicity shown in Eqs.(16) and 
(17) correspond to our particular initial state but it follows essentially from the hamiltonian 
and the cyclic relation (2) and therefore this periodicity is also valid for arbitrary states and 
probabilities. 

From the symmetry of the lattice and of the initial state we expect that the particle will 
diffuse with equal probability clockwise or anti-clockwise, that is, \cn- x {T)\ = \c x (T)\ but 
for the amplitude we may have a different phase on both sides of the initial position. We 
will now show that the amplitude on lattice points symmetric with respect to the initial 
position x = are related by 

c N . x (T) = u~ 2ax c x {T) . (18) 

In order to prove this, consider 

1 N-l 

c N _ x (T) = - E J Ni - p - ]+a) ~ x[p - 1+a) - {p ~ 3)2T ) . (19) 
Using (2) we eliminate N(p — j + a) in the exponent and we add Nx. Therefore 



5 



i N-l 

CN „ X (T) = j-Y: J*v-»-j-°)-<rjrT) . (20) 

p=0 

Now we define another summation index q = N — p with values in {N, N — 1, • • • , 1}. 

i N 

c n - x (t) = j-Y: u (*w-°)-i»-<->rr) . (21) 

9=1 

Since N = 2j + 1, the squared parenthesis in the exponent becomes (q — 1 — j) 2 . Then 

, ,x(l+2j-2oi) JV 

^-,(T) = " Ar £ a;^ 1 -^)-^ 1 "^) . (22) 

Finally, redefining the summation index p = q — 1 and using again (2), we find that the 
right hand side of this equation is u~ 2ax c x (T). A remarkable consequence of relation (18) 
is that, for even N , a quantum particle in a localized state will never diffuse to the antipode 
location. The antipode location, x = N/2, exists only for even N. The proof follows from 
Eq.(18) since we have cn-n/2 = ^~ n ^ 2 cn/2, but u~ N / 2 = —1 therefore cn/2 = —cn/2- That 
is, 

c N /2(T) = VT . (23) 

This is a remarkable result that can be checked by explicit evaluation from Eq.(13) redefining 
the summation index q = p — j running from — j to j. Doing this we obtain a sum whose 
terms are anti-symmetric under q — > — q; therefore they add to zero. Another physically 
appealing proof of this result is provided by Feynman's "sum over paths" method. 4 In this 
case, a path contributing to the probability amplitude for the transition from x = at £j to 
x = N/2 at tf is defined by a set {x^, £&} for each partition of the time interval tj < t k < tf. 
It turns out that for each path {xu, tk} going from x = to x = N/2 there is another path 
{zk\ tk}, symmetric with respect to x — 0, that is, Zk = N — Xk but with the same values of 
{tk}, that cancels its contribution to the probability amplitude, simply because dzk = —dxk- 
In the case N = 2, besides the initial location, there is only one remaining location, the 
antipode. Therefore, for all time, the particle will remain in its initial position. Clearly, for 
N = 2 the states ipo and </?i are not only position eigenvectors but also eigenstates of the 
hamiltonian and therefore they are stationary states. For odd values of N the antipode does 
not exists but we can study the transition probability to diffuse to the "farthest" locations 
x = (N ± l)/2. We will later see a remarkable difference in the odd-iV case. We will see 
that, contrary to what happens in the even case in which the antipodes are never reached, 
if N is odd a sharp distribution will build up in an environment of the antipode at the time 
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T = N/2. This is precisely the time when the state is reconstructed in the iV-even case but 
at the original site. 

We will now calculate the diffusion time, that is, the time that is required for a particle, 
initially localized in one lattice site, to "diffuse" to the whole cyclic lattice. Since the state 
is periodic in time, with period proportional to N, we expect that the state reconstruction 
happens after the whole lattice is visited and therefore the diffusion time should be, at most, 
proportional to N. In order to calculate the diffusion time we must find the time dependence 
of the width of the probability distribution of position. It turns out that for finite N, or for 
periodic distributions, the quantities (^(t),X^(t)} and (*(*), X 2 ^(t)) - (*(*), X^(t)} 2 are 
not appropriate estimates for the center X and width A of the distribution along a cyclic 
lattice or ring. The main reason why they are not appropriate is that any physical quantity 
in a cyclic lattice should be periodic, that is, invariant under x — > x + Na and clearly the 
quantity (\I/(t), X\I/(t)) does not complies to this. Two simple examples: first, let us suppose 
a distribution given by |cat__i| 2 = | Co 1 2 — | ci | 2 = 1/3. Clearly the center of the distribution is 
at the location corresponding to the label x — 0, that is, at the position —aj = —a(N — l)/2, 
but the quantity (${t), X^(t)) is J2x=o a ( x ~ j)\ c x\ 2 — —a(N — 3)/6. For another example, 
consider a uniform distribution that fills completely a ring. In our case of a cyclic lattice 
we have \c x \ 2 = 1/N, Vrr. Clearly, this distribution does not have a center; it should be 
undefined on the ring because all points are equivalent, but the quantity ($(t), X^(t)) is 
EZ£ a(x - j)\c.\* = 0. 

The problem of defining the center X and width A of a distribution in a ring or cyclic 
lattice has been solved 5 ' 6 using the concept of the centroid of a distribution on a ring. Let 
us build a map of the ring into a unit circle in the complex plane. In order to define the 
centroid Z for a probability distribution \c x \ 2 on the sites x = 0, 1, • • • , TV — lofa cyclic 
lattice, let us consider the unit circle in the complex plane with N points located at oo x . The 
centroid of the distribution is a complex number Z = pe %e given by 

JV-l 

Z = pe ie = ]>>1 c *! 2 . (24) 

x=0 

The radial projection of the centroid on the unit circle maps the center of the distribution 
on the ring. Therefore, 

X = a(^N-j), (25) 
and the width A of the distribution is given by 

A 2 = (aN) 2 (l — \Z\ 2 ) , (26) 
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where the factor aN has been chosen such that for a uniform distribution covering the 
whole lattice (Z = 0) the width of the distribution is equal to the size of the lattice. Notice 
that when only one site is occupied, the width is zero. These definitions are clarified in an 
example shown in Figure 2. 

We can now study the time dependence of the width for our initial condition of a particle 
at rest, localized in x = 0. Let us calculate first the centroid. From Eqs.(24) and (15) we 
obtain 

i JV-lJV-liV-l 

Z = — - V ^2 u x u { - p ~ q)x ~ {jp ~ q){ - p+q ~ 2j)T . (27) 
N 2 p=0 q=0 x=0 

The sum over x can be performed: 

JV-l 

E u {p - q+l)x = N(5 g , p+1 + 5 9 , <W-i) • (28) 

x=0 

The first term of the parenthesis corresponds to the vanishing of the exponent of uj and the 
second term is for the case when the exponent is equal to Nx. With the Kronecker deltas 
we perform the sum over q, and the remaining sum over p has a known result. We get then, 



N 



sin 



As expected, the centroid has the same time periodicity as the probability distribution, that 
is, N for odd number of lattice sites and N/2 for an even number of sites. Due to the 
initial condition of a particle in the site x = and to the symmetric diffusion, the centroid 
is real at all times. The study of the time dependence of the centroid, shown in Figs. 3a 
and 3b for N = 16 and 17, allows a simple qualitative description of the time evolution of 
the distribution. In the figures we notice that the centroid oscillates most of the time with 
values close to zero, corresponding to distributions close to (but not necessarily equal to) 
uniform distributions covering the whole lattice. At time T = N/2 the centroid assumes 
the value of Z — 1 in the iV-even case, as expected because at this time the initial state 
is reconstructed, and for odd N it takes the value — (N — 2)/N, close to Z = — 1 for large 
N, implying that at the time T = N/2 the distribution is concentrated at the antipodes of 
the initial location; however, precisely at the antipode there is no lattice site for odd iV and 
the state can not be reconstructed in one location. We see here a sharp distinction in the 
behaviour of diffusion in the even and odd case: the antipode is never reached in the iV-even 
case but the distribution peaks in the neighbourhood of the antipode (at time T = N/2) in 
the iV-odd case. 
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With the knowledge of the centroid, we can now calculate the time dependence of the 
width of the distribution. In particular we want to find the diffusion time T D , that we define 
as the time when the width assumes its maximal value aN for the first time. Notice that 
when the centroid vanishes, the width assumes its maximal value. From Eq.(29) we see that 
the centroid vanishes for T = 1 for all N, therefore the width is maximal (aN) at T = 1. 
However Eq.(29) has another root for a time T smaller than 1 when N > 4. Of course, when 
N = 2 the diffusion time is infinite because the particle never diffuses out of the initial site. 
Summarizing we have 

{oo , for iV = 2 
1 , for N = 3 . (30) 
K for TV > 4 
2(JV-2) ' 1UI iv ^ ^ 

It might at first seem strange that the defined diffusion time decreases towards a constant 
value T D = 1/2 with an increasing number of sites N, but we can see that this is to be 
expected as a consequence of indeterminacy principle. Increasing the number of sites iV with 
the same initial condition of a particle in one site, is equivalent to a sharper localization of the 
initial state. This implies a wider momentum spread, responsible for a faster diffusion that 
decreases the diffusion time. The explicit time dependence of the width of the distribution 
is then given by 



{N \> } +1) • (31) 




~ \| V sin (f T ) 

This quantity is zero at T = 0, grows with time, and takes the maximal value aN at T = To; 
then it oscillates with values close to the maximal value except at time T = N/2 when the 
width becomes zero for iV-even or decreases to a2y/N — 1 for odd N. At this time, that we 
call first reconstruction time Tr = N/2, the particle is reconstructed at the original site (N 
even) or is concentrated near the antipode (N odd). Notice that at this reconstruction time 
Tr the state is reconstructed only if N is even whereas for odd N the probability distribution 
for the location of the particle peaks but there is no exact reconstruction of the particle in 
one location of the antipode. For very short times T <C Tp the system does not notices the 
geometry of the cyclic lattice and the width grows linearly with time with a diffusion speed 
increasing with the lattice size N. Indeed, the first term in the Taylor expansion of A is 



A = a2n^\(N - 1)(N - 2) T for T < T D . (32) 

We can now investigate whether the reconstruction of a localized state for the particle 
at time T R = N/2 at the original site (N even) or the concentration of the particle near the 
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antipodes (N odd) is affected by the parity of the initial state. The initial state considered 
above, a particle in one site, has necessarily even parity. In order to be able to study also 
the effect of an odd parity initial state we will consider an initial state of a particle at 
rest, (P) = 0, in an even or odd superposition of two neighbouring position eigenstates 
corresponding to the sites x = and x — 1: 

*±(0) = ^ (voiw'Vi) (33) 

With this initial state, we can calculate the time evolution of the centroid. However, the 
centroid will no longer be a real number. It is therefore convenient to make a rotation of 
the centroid in the complex plane by an angle co- 1 / 2 in order to obtain the real quantity 
Z±(T) = uo~ l l 2 Z±(T), where Z±(T) is the centroid corresponding to the two initial states 
*±(0). This results in 

l( ^ sin (j^V - 1)2T) sin (%(N - 1)(2T - 1 
Z±{T) = — cos — , 7- — x ± 



iV 



sin(^2T) 2sin(f(2T-l 



sin(f(7V-l)(2T + l)) /tt. . 



2sin(f (2T + 1)) ' \N, 

In Figures 4a and 4b we see the time evolution of the (rotated) centroid Z + {T) for an even 
initial state ^+(0) for even and odd N. For a qualitative comparison with Figure 1, we have 
taken N = 33 and N = 34 in order to have similar relation between the size of the lattice and 
the number of sites of the initial state. From this comparison it is clear that the behaviour 
is similar. At time T = N/2, a localized even state is reconstructed at the original locations 
if iV is even or the particle is localized at the antipodes if iV is odd. In Figures 5a and 5b 
we can see that this is also true when the initial state \P_(0) is odd, but the effect is much 
blurred by rapid oscillations of the centroid. Shortly before and after every reconstruction 
of the particle, it is almost reconstructed but on the opposite side of the lattice. For both, 
even and odd parity states, the initial value of the centroid, Z±(0) = cos(n/N), is exactly 
recovered for even N at time T = N/2 (this must be so because the state is periodic) and 
for odd N the centroid reaches the minimum value Z±{N/2) = — ((N — 2) cos(7r/iV) ±2)/N. 
For large N this minimum value approaches — cos(7r/iV), corresponding to the occupation 
of two neighbouring sites at the antipodes. 

IV. THE CONTINUOUS LIMIT 

We have found that there are very strong differences in the behaviour of the system when 
iV takes even or odd values. Of course, all these differences must be compatible with the 
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continuous limit when iV — > oo where we can not differentiate between even or odd N. In 
this section we will investigate this limit. First we must redefine the indices of summation 
in a symmetric way such that they can take positive and negative values. Let then 

y = a(x-j) G [-aj,aj] 

q = g(p-j) e [-gj,gj] ■ 

Anticipating that in the limit A" — > oo the position and momentum eigenfunctions will not 
be normalizable, we define these eigenfunctions in terms of the symmetric indices as 

1 1 

(f y = —= (p x and (p q = — (j) p . (36) 

V a V9 

If in the limit A" — > oo we also take a — > or g — > then the summations become integrals 
according to the scheme 

a J poo 9 J poo 

J^a-^ / dy or -> / dg . (37) 

The limit A" — > oo is constrained by the condition Nag = 2n and therefore we will 
consider three different limits LI, L2, L3, that will correspond to three different physical 
systems: 

f N -> oo, a -y 0, o -> 0, 
LI : < ' ' y ' (38) 

lye [— oo, oo], g G [— oo, oo] . 



N -> oo, a -> 0, Na = L, g = \ 
ye[-L/2,L/2], q=^n, n 




L2 -{ r r ,o r 2. [ ±V2, ±3/2, • • • iVeven, (39) 

■ ■ Nodd. 



N -> oo, y -> 0, Ng = G, a = §, 
L3 : 4 ^ ^ 2^ f i 1 / 2 , ±3/2, • • • A^even, (40) 



?6 [-67/2,(7/2], i/ = fn, 



0,±1,±2,--- Nodd. 



In the limit LI both variables y and g are continuous and unbound whereas in L2 the variable 
y is bounded and continuous but q is unbound and discrete; these properties are exchanged 
in L3. 

In the limit LI, the physical system becomes a free particle moving in a one dimensional 
infinite space where position and momentum observable can take continuous values. In the 
limit L2, the physical system is a free particle moving in a ring of perimeter L. Position is 
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continuous and takes values from — L/2 to L/2 whereas momentum is a discrete variable. 
We will later see that among the two choices for the number n, only the values 0, ±1, ±2, • • • 
are physically meaningful. This system also corresponds to a particle in a box with periodic 
boundary conditions. Finally, in the limit L3, the physical system is a particle moving in a 
one dimensional infinite lattice with lattice constant a = 2-k/G and continuous momentum 
restricted to the Brillouin zone [-G/2, G/2\. 

The striking differences in the behaviour of the system between even and odd N appear 
in the time periodicity of the state and probability and in the first reconstruction time for 
the probability distribution. These differences involve a time scale t = Nr oc Na/g = 2n/g 2 . 
In both limits, LI and L3, this time scale is infinite and therefore we should not worry about 
whether N is even or odd when taking the limit iV — > oo; however in the limit L2 the time 
scale is finite and proportional to L 2 . In this last case we will see that the even N case is 
mathematically sound but does not corresponds to any reasonable physical system. 

It is convenient, in order to analyse the LI and L2 limits, to adopt the position repre- 
sentation of the eigenfunctions where the momentum eigenvectors are given by Eqs.(7,36) 
as 

Uv) = fa, M = -L e^+^y-^ . (4i) 

v Zn 



In the limit LI, where a — > and g — > 0, this eigenfunction becomes 

My) = ^= ^ , (42) 

provided that, in the even N case (a = 1/2), The values of y and q remain finite (otherwise 
a minus sign can appear). Expanding the position eigenfunctions in the momentum basis 
we obtain 

My) = ^ f°°dq ^ y - yl) = S(y - y') . (43) 

We obtain therefore the usual position and momentum eigenfunctions for a free particle 
moving in a line. 

Let us now consider the L2 limit where we have two possibilities: a = 0, n = 
0, ±1, ±2, • • • and a = 1/2, n = ±1/2, ±3/2, • • •. In the first case Eq.(41) results in 

<t> q {y) = -L e^ n , n = 0, ±1, ±2, • • • , (44) 
v Zn 



and in the second case, assuming m = ±1/2, ±3/2, • • •, we have 

= _L e iy 2 -tm+i^y = _L e iv£(m+§) = _L e *V^n ? n = , ±1 , ±2, • • • . (45) 

\/2tt \/27t \/2ti 
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Therefore both cases lead to the same position representation of the momentum eigenfunc- 
tion. So far it would seem that the even and odd N cases are identical in the limit N — > oo 
however this is not so as we will see next. It turns out that in the even N case, when 
a — 1/2, the momentum operator in the position representation is not given by the deriva- 
tive operator as usual. In order to prove this, consider the first equation in (10) written in 
terms of the symmetric variables, that is, 



e-** <p y = e^ a V[y+a] = e^ aa <p [y+a] . (46) 
Applying the limit L2 we get, 

(1 - iaP) tp y = (1 + i—aa) (p y+a . (47) 

Therefore 

P Vy = t lim — - -j-a ip y+a , (48) 

where we see that only in the odd iV case, when a = 0, is the momentum operator given by 
the derivative operator. 

The inadequacy of even iV in the limit is more conveniently seen if we absorb the phase 
e ia(gy-aq) j n eigenfunctions as was mentioned at the end of section II. In this case the 
a-dependent phase in Eq.(41) would not appear and we would have two different position 
representations of the momentum eigenfunctions given by 

<f>Uy) = -L e iw T(o. ±L ±2 > -) for odd N 

tffy) = J— e W 2 -?(±l/2, ±3/2, -..) fo even N 



The momentum eigenfunctions 4>\{y) are the same as the ones obtained before in Eq.(44) 
and the other ones, </>g(y), are mathematically sound but are inadequate for physical systems 
because they are anti-symmetric, 4> 2 q {— L/2) = —(p 2 AL/2) and have period 2L, whereas all 
reasonable physical states for a particle in a ring are symmetric and have space periodicity 
L. 

As a further confirmation that the L2 limit corresponds with the odd case, we will 
show that an initial state in a ring is reconstructed at the antipodes at the reconstruction 
time tn = T^t = Nt/2 = mL 2 /(2ir) as it happens in the case of finite but odd N . In order 
to prove this we assume an arbitrary initial state expanded in terms of the momentum base 

ij(y,0)=Y,c Q My) ■ (49) 
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We apply the time evolution operator to this state, considering that 



and using Eq.(44) we get 



M,*) = 4= E c n e-^(^) V ' + ^ . (51) 

V 



n=0,±l,±2,- 

2 



Consider now this state at the reconstruction time tR = mL /(2n), 

^y,t R ) = -^= J2 c n e-*™ 2 e^ n . (52) 



2lT n=0 ,±l,±2 / 



Now, since n 2 and n have the same parity, it is e mn = e mn , and we get, 

1 

27 



lKl/,f*) = 4= £ c n e^)^. (53) 



n=0,±l,±2,- 

Therefore 

^(j/,**)=^(j/-^,0) . (54) 

with the meaning that the state at the reconstruction time is equal to the initial state, 
but shifted to the antipode y — L/2. 

A particle in a localized state in a continuous ring will flip back and forth between the 
original position and its antipode. If the particle is electrically charged, the system will 
radiate electromagnetic energy and decay to a nonlocalized stationary state (not necessarily 
the ground state) and a small fluctuation away of a nonlocalized state will grow, absorbing 
electromagnetic energy of the appropriate frequency (for an electron in a ring of 10/i to 100/i 
of perimeter, the radiation will be in the radio frequencies). This radiation, or the corre- 
sponding absorption, could be experimentally detected, in particular if one builds a material 
with a large number of conducting rings. Such a material, whose dielectric properties follow 
from a fundamental quantum mechanical effect, could find technological applications. 

Finally, the L3 limit is treated equal to the L2 case but in terms of the momentum 
representation of the eigenfunctions. Similar arguments show that the even N case leads, in 
the limit, to unphysical situations. 

V. CONCLUSION 

In this work we have studied the diffusion of a quantum mechanical particle, initially 
localized, in a ring with iV sites. This diffusion has qualitative features quite different from 
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the diffusion of a particle performing a classical random walk. It is well known that in 
a classical random walk, the width of the distribution grows like y/T whereas quantum 
mechanical diffusion grows initially proportional to T. Furthermore, we see in Eq. (32) that 
the speed of quantum diffusion, for large N, increases linearly with the size of the lattice 
aN. This non local effect is contrary to the classical behaviour and can be understood 
qualitatively as a consequence of Heisenberg's indeterminacy principle: if the initial state 
is a particle in one site of the lattice, increasing the number of sites is equivalent to a 
sharper localization relative to the lattice size, and this results in a wider momentum spread, 
responsible for the increase in diffusion speed. Since the diffusion speed increases with the 
number of sites N it is reasonable to expect that the time necessary to diffuse to the whole 
lattice will be constant independent of the lattice size. This is indeed the result shown in 
Eq.(30) where we see that, for large N, the diffusion time Tp is constant. This is again in 
contradiction with the behaviour of the classical random walk where the covering time 7 (the 
time it takes for a random walk to visit all the lattice sites) for a cyclic lattice increases 
quadratically with N (precisely, N(N — l)/2). 
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FIGURE CAPTIONS 



FIGURE I. Cyclic lattice with N sites characterized by a label x running from x = to 
x = N — 1 and lattice constant a. The position observable corresponding to site x has the 
eigenvalue a(x — j) where j — (N — l)/2 and can take positive and negative values. 

FIGURE 2. An example for the centroid of a distribution. The symbol © shows the position 
of the centroid Z = pe ld for a distribution where the filled dots have a constant occupation 
probability and all other sites are empty. The center X of the distribution is shown and the 
width A is proportional to the chord C shown. 

FIGURE 3. Time dependence of the centroid for even (a) and odd (b) number of sites. At 
time N/2 the state is reconstructed at the original site for even N and is concentrated at 
the antipodes for odd N . 

FIGURE 4. Time dependence of the (rotated) centroid for even (a) and odd (b) number of 
sites for an initial even state occupying two neighbouring sites. At time N/2 the state is 
reconstructed at the original sites for even and is concentrated at the antipodes for odd 
N. 

FIGURE 5. Time dependence of the (rotated) centroid for even (a) and odd (b) number 
of sites for an initial odd state occupying two neighbouring sites. At time N/2 the state is 
reconstructed at the original sites for even N and is concentrated at the antipodes for odd 
N. 
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